Continuum Mechanics and Thermodynamics vol. 24 pp. 505-513 (2012). 



H. Gouin • A. Muracchini • T. Ruggeri* 



(N 

o 

(N 

> 
O 



Of 



> 
o 

(N 



(N 



X 



On the Miiller paradox for thermal-incompressible 
media 



Abstract In his monograph Thermodynamics, I. Miiller proves that for incompressible media the 
volume does not change with the temperature. This Miiller paradox yields an incompatibility between 
experimental evidence and the entropy principle. This result has generated much debate within the 
mathematical and thermodynamical communities as to the basis of Boussinesq approximation in fluid 
Oi dynamics. 

\ The aim of this paper is to prove that for an appropriate definition of incompressibility, as a limiting 

case of quasi thermal- incompressible body, the entropy principle holds for pressures smaller than a 
critical pressure value. The main consequence of our result is the physically obvious one, that for very 
large pressures, no body can be perfectly incompressible. The result is first established in the fluid 
c/2 " case. In the case of hyperelastic media subject to large deformations the approach is similar, but with 

O , a suitable definition of the pressure associated with convenient stress tensor decomposition. 

Keywords Incompressible fluids and solids • Entropy principle for incompressible materials • 
^ <—j Boussinesq approximation 



PACS 47.10. ab • 62.20.D- • 62.50.-p ■ 64.70.qd ■ 65.40.De 

Mathematics Subject Classification (2000) 74F05 • 74A10 • 74A15 • 76Bxx • 76E19 
1 Introduction 

It is well known that compressible and incompressible bodies have different mathematical treatments: 
for compressible media the pressure is a constitutive function, while for incompressible media the 
pressure comes from a Lagrange multiplier associated with the constraint of incompressibility. From 
an experimental point of view, incompressible medium has no real existence but can be approximated 
as limit case of the compressible one. Starting from this observation much literature has been devoted 
by using qualitative analysis and numerical methods to search solutions of the incompressible case; 
for example, the limit of solutions of the compressible fluids are considered as the Mach number tends 
to zero under certain assumptions on the initial data (see e.g. the isothermal case [1,2,3,4 ). When 



the thermal effects are relevant the limit is more ambiguous and depends on the particular model of 
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incompressibility (see e.g. [5]). With the aim to have the same set of equations for compressible and 
incompressible fluids, Miiller [5J and others (see e.g. [TIE]) choose the pressure p as unknown field 
variable instead of the density p. More precisely, in the case of dissipative Navier-Stokes-Fourier fluids, 
they add to the balance law system of mass, momentum and energy, 

dp , dm = Q) (i) 



dt dx 
dpvj i d 



(2) 



dt + ~dx~y ViV] ~ tlj) = ^ 

(4) 

d P {e+\v 2 ) d f / 1 2 \ 1 

p £ + -« ] Vi - UjVj + qi > = pfjVj + pr, (5) 



A/ 



at 9a;, r V 2 

the following constitutive equations, 

p = p(p,T), e = e(p,T), (6) 
<7ij = ^(p, T) dfe fe % + 2/x(p, T) dy, g, = -n(p, T) dT/dxi , (7) 

with 

t = -pi + a , (8) 

where t = (ty) represents the stress tensor, rr = (cry) is the viscous shear stress, 1 = (5y ) the identity 
tensor and Sij the Kronecker symbol. The other variables have the usual meaning: v = (vj), q = (qj), 
e, T respectively denote the velocity, the heat flux, the internal energy and the temperature while 
f = (/ 3 ), r, are the specific body force and the heat supply; k is the heat conductivity and the scalars 
v and p are the viscosity coefficients. Matrix d = \\dij\\ denotes the symmetric part of Vv (dkk = div v). 
Consequently, system ©-© is a closed system for the unknown variables p, T, v. 
Miiller defines an incompressible fluid as a medium for which the constitutive equations (O-© are 
independent of the pressure, and in particular: 

p = p(T), e = e(T). (9) 

Nonetheless, he proves [B] that the only function p(T) compatible with the entropy principle is a 
constant function p — p a . Obviously, this result disagrees with experiments showing that the density 
changes with the temperature (see e.g. [5]) and with all the theoretical results given in particular in the 
so-called Boussinesq approximation (see e.g. QUE])- We call this contradiction the Miiller paradox. 
In this paper we prove that for a convenient definition of incompressibility, the entropy principle is 
compatible with the fact that volume changes with temperature provided that the pressure is smaller 
than a critical pressure value p cr . These results provide a precise value of the critical pressure under 
which a fluid can be experimentally similar to an incompressible one and permits to obtain a quanti- 
tative measurement for which the Boussinesq approximation can be considered as valid. 
The paper is organized as follows: 

In Section 2, the Miiller paradox is examined in fluid case. 

In Section 3, we prove that exists a critical pressure p cr such that for any p <C p cr the paradox is 
solved. Then, a numerical evaluation is done allowing to obtain the critical pressure value for which 
water can be considered as incompressible liquid. 

In section 4, the case of hyperelastic media subject to large deformations is considered; thanks to a 
stress tensor decomposition carrying out a suitable definition of the pressure, an approach similar to 
fluid case allows to solve the Miiller paradox. We perform a numerical evaluation of the critical pressure 
for pure gum rubber. 
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2 The Miiller paradox for incompressible fluids 

The Gibbs relation 

TdS = de Trdp, 

P 2 

where S is the entropy density, can be rewritten by using the chemical potential 

fj, = s+pV -TS, (10) 

in the form 

dp = Vdp- SdT , (11) 

where V — 1/ p is the specific volume. The choice of the chemical potential is natural if we use p and 
T as variables. In fact from eqs. (|10till[) it follows 

V = fi p , S = -a T , e = a - pp p - Tp T , (12) 

where 

dp\ ( dp 



Mp= [dp-J T ' ^ \dT 
From eqs. (fHJh - 3 we get 

£ P = -PV P - TV T . (13) 
The Miiller definition of incompressibility ([5]) and Eq. (fT5T) imply V is constant, i.e. 

p = Pa = constant . 

Obviously, this Miiller result [5] (p. 27) disagrees with the experiments proving that the volume of 
fluids changes with the temperature and also, a little, with the pressure [51IT2"]. 

3 Removal of the Miiller paradox for fluids 

According to the fact that the so-called compressibility coefficient j3 = —V p /V is very small in the 
case of incompressible body, to remove the Miiller paradox we define the notion of quasi thermal- 
incompressible fluid: 

A quasi thermal-incompressible fluid is a medium for which the only equation independent of p 
among constitutive equations (O-© is the density: 

P = P(T) • (14) 

We assume that the condition of independence of p is not necessary for the other constitutive equations; 
in particular, the internal energy e remains function of p and T: 

E = e(p,T). 

From Eq. (|14p . by integration of Eq. (fT2"|) i . we obtain 

n(p,T) = V(T)p + iio(T), (15) 

and substituting in Eq. (TT2"lU 



e(p,T) = -TV'(T)p + e(T), (16) 
where 

e(T) = (j, - ?>o , 

with ' = d/dT. Following Miiller proposal, we name incompressible fluid (or better perfectly incom- 
pressible fluid) a fluid for which all the constitutive equations are independent of p. 
Therefore, from Eq. (I16[) , a quasi thermal-incompressible fluid tends to be perfectly incompressible if 
e(p,T) can be approximated with e(T), i.e. when: 

e(T) p^T) 

P<< WW = 1W (17) 
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For example, let us consider the classical linear behavior (also typical of Boussinesq approximation): 

V = V [l + a(T-T )], (18) 

where a is the thermal expansion, the constant Vb = 1/po is associated with the scale of volume and Tq 
is a reference temperature. If we assume that the specific heat at constant pressure c p = hx is constant 
(where h — e + pV is the enthalpy), then from Eq. (fTl)]) we have 

e{T) = c p T. (19) 

Inserting eqs. ()18til9|l into inequality (fTT)) we obtain: 

P<.Pcr with p cr — c p po/a. (20) 

We call p cr the critical pressure at density po- As we expect in physical situations, we note that p cr 
is inversely proportional to a. Critical pressure p cr and inequality (|20[) characterize the fact that a 
quasi thermal-incompressible fluid is experimentally similar to a perfectly incompressible fluid. In such 
a case, the Miiller paradox is removed. 



3.1 Application to water 

To evaluate the magnitude order of the critical pressure (|20p . we give some numerical results in the 
case of water. 

At temperature T = 20° C, we get [51IT5] : 

po - 10 3 kg/m 3 , c p ~ 4.2 ■ 10 3 Joule/kg.K, a = 207 • 1(T 6 / K > 
and from Eq. (20]) we deduce: p cr ~ 2 • 10 10 Pascal ~ 2 • 10 5 atm. 

The value of critical pressure is large with respect to the normal pressure conditions. The fact that for 
usual pressures, a liquid is experimentally incompressible and the volume changes with the temperature 
( Eq. (|18|) ) does not violate the principles of thermodynamics. We can interpret the result in another 
way: for very large pressure, perfectly incompressible fluids do not exist. From physical point of view, 
this observation seems reasonable, for example in astrophysics where very high pressures are present. 
We also observe that from experimental data of sound velocity in the water [121113] , the compressibility 
coefficient (3 is not zero but very small (j3 — 4.98 x 10 -10 / Pascal). This is in agreement with the fact 
that perfect incompressibility is an idealization. 

Remark 

Our definition of quasi thermal-incompressibility is, for some authors (e.g. (7]|5]), the definition of 
incompressible body in the context of Boussinesq approximation. Nevertheless, we strongly believe 
that the Miiller definition of incompressible fluid is the correct one even if it is a limit case: only in this 
case, we obtain the same differential equations as for the usual incompressible approach with p, v, T 
variables. In fact in the case of incompressibility, i.e. when p = p(T), a generic constitutive quantity 
4> = <f>(p, T) becomes a function only of T. In our analysis perfectly incompressible fluid is considered 
as limit case of quasi thermal-incompressible fluid. 

The quasi thermal-incompressibility is obtained as a limit process justifying the compatibility between 
incompressibility and Gibbs relation when inequality (|20|) is verified. The quasi thermal-incompressibility 
does not characterize a real compressible material; for real compressible fluids the chemical potential 
p must be a concave function of (p, T). When V depends only on T, the chemical potential is a linear 
function of p (see Eq. (|15p ) and consequently cannot be concave. For real compressible materials the 
volume V necessarily depends on p; quasi thermal-incompressible materials can be considered as an 
approximation of incompressible materials when the pressure is sufficiently small such that inequality 
is satisfied. 

Moreover, note that we use term of quasi-thermal incompressibility, rather than term of quasi incom- 
pressibility, because this last term has a different meaning in the pure isothermal mechanical case. 
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4 Incompressible hyperelastic media 

In the case of elasticity, Miiller presented a similar paradox [B] (p. 263). In fact, he proved that 

J = detF, 

where F denotes the deformation gradient, cannot depend on the temperature T and must be constant 
as in the pure mechanical case, i.e. J = 1. This result also disagrees with experiments proving that the 
volume of an elastic incompressible solid changes with the temperature. 

We use similar arguments as in fluid case, but on the contrary of the Miiller procedure, we consider a 
particular decomposition of the stress tensor allowing to remove the paradox. Moreover, our approach 
is also valid for non-isotropic materials. 



4.1 A decomposition of the stress tensor in the case of hyperelastic medium 

Each particle of the continuous medium is labeled by a material variable X, ranging from a reference 
configuration Dq into an Euclidian space |14j . The reference density po is given as a function on T>q 

m 

The expression x = </>(X, t) of the spatial position describes the motion of the continuous medium. 
Generally, </>(., t) is a twice continuously diffcrentiablc diffcomorphism from T>q into a compact oriented 
manifold T> t constituting the image of the material at time t. As usual, we denote by C = F T F the 
right Cauchy-Green deformation tensor, where superscript T means the transposition. Let us recall: 

p (detC)^ =p - (21) 

For the sake of simplicity, we consider the case of homogeneous bodies (the case of inhomogeneous 
bodies can be treated in the same way) . The internal energy density is supposed to be function of the 
tensor C and temperature T 

e = s{C,T), 
and the stress tensor can be written [TMTMTT] . 



2 P f|^F t , (22) 



where ip = e — TS is the specific free energy. By writing [TMlU] : 



C=(^) C, (23) 



(detC)3 V P 

the specific free energy can be expressed in the form: 

tl> = f(p,C,T). 

We note that C is substituted by the independent variables p and C. Since detC = 1, the variable 
p corresponds to the change of volume while the tensorial variable C represents the distortion of the 
medium. This point is fundamental for the decomposition of the stress tensor and will be the key of 
the demonstration. When / is independent of C we are back to the fluid case. It is more convenient 
to introduce the function g such that: 



-'"■ c ' T)s/ l"'(3^c)I C ' r 

Consequently, g is a homogeneous function of degree zero with respect to C. 
From Eqs. (|2HI22p it follows: 
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Differentiating Eq. (12T1) and using Jacobi's identity, we obtain 

dp = -\pC~ l -dC, 

where the dot represents the scalar product between matrices. Hence, 

dp 1 



and consequently, Eq. (jM)) yields, 



2 d 9 , , o^ d 9^T 



t = - p f P 1+2pF ob F 

Due to the fact g is homogeneous of degree zero with respect to C, from the Euler's identity we 
immediately deduce: 

c-o 

and 

t = -pl+T, (25) 

with 

P=p2 % i T = 2pT §h FT and trr = ' ( 26 ) 

where tr is the trace operator. Let us note that, in relation (|25p. t is similar to the fluid decomposition 
of Eq. flSJ and in the solid case, p is analog to a pressure. 

The decomposition (|25H26[) allows to define a pressure also in the case of an elastic body. 



4.2 Removal of the Muller paradox for hyperelastic media 
The Gibbs equation in the case of elastic materials is P~51fT^lfT 7 7) 

TdS = de- — S • dC . (27) 
2po 

where 

S = JF~ 1 t(F T y 1 (28) 

is the second Piola-Kirchhoff stress tensor. Inserting Eq. (j25|) into Eq. ((28]) . we get from Eq. (|27|) the 
Gibbs relation: 

TdS = de-^dp--^-F- 1 rF- lT ■ dC. (29) 
P 2p 



From Eq. (|23l) we obtain: 



d C = -~- dp+(^\ dC, 

3 p V p 



and Eq. (gSJ) leads to: 

2 

TdS = de-^dp+ -i^F-V F- 1 " • C dp - — ( ^) F-V F- 1 " • dC. 
P 2 3 p 2 2 p V P ) 

From Eq. (|2SJ) , we get F-VF- lT • C = tr (f-^F-^c) = trr = 0. If we define 
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we obtain the Gibbs relation in the final form: 

TdS = de - 4 dp + r ■ dC. (31) 
P 2 

As in the case of fluid, introducing the chemical potential 

p = e + pV + r-C-TS, 

Eq. (JHU) implies 

dp = Vdp - SdT + C-df. (32) 

Equation (|3T)f together with Eq. ([2^)1 imply tr(rC) = 0. Therefore r has only five independent 
components. 

From Eq. (|3"2")l it is natural to introduce the change of variables (p,T,C) into (p, T, t). In fact, Eq. 
((32)1 implies: 

V = p p , C = p? , S = -p T , £ = p - Tp T -pp p -T ■ pt , (33) 

where now 

dp\ ( dp\ ( dp~ 



By analogy with the case of fluids, a quasi thermal-incompressible elastic medium must verify 
the condition: 

J = J(T) or equivalently V = V(T) , 
while the other constitutive equations still depend on (p,T, t). In particular 

e = e(p,T,T). 

By integration of Eq. ([33| i , we obtain: 

p{p,T,T) = V(T)p + po(T,T), 
and by substituting in Eq. (j33|) zi. we get 

e = -TV'{T)p + e(T,T), (34) 

with 

e(T, f) = p - Tpor - t ■ p Q f . 

Following Miiller proposal again, we define an incompressible elastic medium (or better perfectly in- 
compressible elastic medium) as a solid for which all the constitutive equations are independent of p. 
Therefore, taking account of Eq. (|34[) . a thermal-incompressible elastic medium tends to be perfectly 
incompressible if s(p, T, r) can be approximated by e(T, r). This assumption is verified when 

e(T, t) 

p<< wiy- 

We consider the case: 

e{T,T)=c p T, V(T) = V [l + a{T-T )}, 

where now, a and Vo are positive constants and c p can be a function of t. As well as for fluids, we 
obtain: 

P < Cppo/a. (35) 

We assume that c p po/a, varying with r, has a minimal value p cr ; we again denote by p cr the critical 
pressure. In this case the critical pressure p cr and inequality p <C Per are characteristic of quasi thermal- 
incompressible elastic medium to be experimentally similar to perfectly incompressible one and the 
Miiller paradox is removed. 
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4.3 Application to pure gum rubber 

The most famous incompressible hyperplastic medium is the pure gum rubber. This material was 
studied by many authors and in particular by P.J. Flory, Nobel prize of Chemistry in 1974 |20j . 
In the range of temperatures [50° C, 85° C] physical constants are [T3"] : 

po ~ 930kg/m 3 , c p ~ 1.9 • 10 3 Joule/kg.K, a = 6.7 • 10 _3 /K. 

Equation ([33)) allows to obtain: p cr ~ 2.7 • 10 s Pascal ~ 2.7 • 10 3 atm. The critical pressure is very 
large and therefore pure gum rubber is a good incompressible body in normal conditions with usual 
pressures. 

5 Conclusions 

In the case of fluids and elastic media, we showed that, for pressures smaller than a critical value, the 
volume can depend on the temperature if incompressibility is defined as limit case of quasi thermal- 
incompressibility. Quasi thermal-incompressible materials can be considered as an approximation of 
incompressible materials - in the sense of [TIE] where j3 — - if the pressure is small enough such that 
inequality (|2"tJ)) is satisfied. 

To obtain these results, we used temperature as natural thermodynamical variable. Nevertheless, Man- 
acorda first noted (see also [2"2"ll23| ) that in the case V = V(T), instabilities occur in wave prop- 
agations. The instabilities are due to the chemical potential non-concavity (see Remark in previous 
section) and the sound velocity c becomes complex. For this reason, some authors consider the volume 
as function of entropy V = V(S) instead of function of temperature p4ll25] . This assumption does not 
seem realistic: V — V(S) cannot be measured because entropy is not an observable and moreover, in 
this case 1/c — (i.e., the sound velocity is infinite); as a consequence the mathematical structure of 
Euler fluids becomes parabolic. 

Our goal was to present the simplest model for removing the thermodynamical paradox when the 
volume depends only on T. Nevertheless, a more realistic definition of quasi thermal-incompressibility 
needs to suppose that the compressibility coefficient f3 is small but not zero; the concavity of the chem- 
ical potential can be restored and the sound velocity can be real with consequence that incompressible 
body can be seen as a limit case of a compressible one but the present result cannot quantitatively 
change. This will be the subject of a forthcoming study. 

Finally thanks to decomposition (|25l) . the technique used to remove the Muller paradox is available 
both for fluids and elastic media. This similarity allows to forecast a possible Boussinesq approximation 
in case of elastic media. 

Acknowledgments: The authors are grateful to Professor Salvatore Rionero for his interest on this paper and 
his useful reference suggestions. 
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